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Abstract. We study the importance of interband effects on the orbital susceptibility of three 
bands a-Ts tight-binding models. The particularity of these models is that the coupling between 
the three energy bands (which is encoded in the wavefunctions properties) can be tuned (by 
a parameter a) without any modification of the energy spectrum. Using the gauge-invariant 
perturbative formalism that we have recently developped [1], we obtain a generic formula of 
the orbital susceptibility of a-Ts tight-binding models. Considering then three characteristic 
examples that exhibit either Dirac, semi-Dirac or quadratic band touching, we show that by 
varying the parameter a and thus the wavefunctions interband couplings, it is possible to drive 
a transition from a diamagnetic to a paramagnetic peak of the orbital susceptibility at the 
band touching. In the presence of a gap separating the dispersive bands, we show that the 
susceptibility inside the gap exhibits a similar dia to paramagnetic transition. 


1. Introduction 

The orbital magnetic susceptibility of free electrons was computed long ago by Landau [2] 
and was found to be diamagnetic. Subsequently, Peierls extended Landau’s result to the case 
of a single band tight-binding model. He found a formula for the orbital susceptibility that 
only depends on the zero-field band energy spectrum. This result already showed that the band 
structure can have a strong influence on the magnetic response of crystals [3]. For example, near 
a saddle point of the dispersion relation, the Peierls orbital susceptibility becomes paramagnetic 
[4]. However, one important effect was left out, namely the coupling between the bands in the 
case of several bands. A striking example is the possibility of having a finite orbital susceptibility 
inside the gap of a band insulator at zero temperature. This was understood as an inter-band 
effect by Fukuyama and Kubo on the example of bismuth [5]. Fukuyama also provided a compact 
and quite general linear-response formula that includes interband effects [6]. However, despite 
its many successes, this formula does not work for non-separable tight-binding models. The aim 
of this paper is to present orbital susceptibility results obtained from an exact linear response 
formula that we recently derived for tight-binding models [1]. In order to show the importance 
of inter-band (or band coupling) effects, we consider here a family of tunable three band tight- 
binding models constructed on a Ts lattice (or dice lattice) and that we call a-Ts [7]. The 
dice model was first introduced [8] as a simple 2-dimensional model that exhibits localized 



and extended states as the same time. It also presents some peculiar properties under strong 
magnetic field [9]. These a-Ta models depend on a real parameter a and have the important 
property that their zero-field energy spectrum -which is essentially that of graphene with an 
additional flat band- is independent of a. However, the parameter a has a strong influence on 
the zero-field eigenstates and therefore on the Berry curvature. This influence is revealed by the 
orbital susceptibility that changes sign as a function of a. 

The paper is organized as follows. In §2 we present the three bands a-Ts tight-binding models 
and characterize the key properties of their energy spectrum and wavefunctions. In §3, using 
the gauge-invariant perturbative formalism that we have recently developped [1], we provide a 
generic formula for the orbital susceptibility of a-Ts tight-binding models. In §4 we apply this 
susceptibity formula to three characteristic examples of a-Ts tight-binding models that exhibit 
respectively Dirac, semi-Dirac and quadratic band touching at low-energy. In §5 we summarize 
the main results of this study. 


2. Tight-binding models on the Tjj lattice 

Starting from the honeycomb lattice with two sites {A,B) per unit cell, the Tjj or dice lattice 
is obtained by connecting additional (C) sites at the center of each hexagon to the B sites (see 
Fig. 1). The dice lattice is thus a triangular Bravais lattice with three sites {A,B,C) per unit 
cell. We consider tight-binding models that consist of spinless electrons hopping on this lattice. 
In its simplest form, we allow for a constant onsite potential term -|-A on sites A, C and —A 
on sites B and an isotropic nearest-neighbors hopping with amplitude Cat from A to B and Sat 
from C to B with Cq = ^ such that -|-s^ = 1. The real space representation 

of the corresponding Hamiltonian follows as 


h — [^a't{SrB-Si,rA ^rB-S2,rA + <^rB-(53,rA) I'^s) | 

~\~Sat{6rg-\-Si,rc T ^rB+S2,rc T ^‘rB+S3,rc^\'f' b) {fcW T h.C, 

where di,52)<53 are the three vectors connecting nearest-neighbors sites. We assume that 
the localized orbital basis is orthogonal such that the position operator is 

purely diagonal (r = s c Iiitroducing the Bloch states basis \kj=A,B,c) = 

YZvj |rj), for each wavevector k = {kx, ky), the Bloch Hamiltonian matrix associated to this 
a-Ts tight-binding model reads 


/ A Cafk 0 \ 

hk = Cafl -A Safk , (2) 

V 0 Safi A / 

where fk = In the following, we will also consider 

generalized versions of this model that essentially consist in a modified fk- This class of models 
interpolates between the honeycomb {a = 0) and the isotropic dice lattice (a = 1). 

The remarkable and interesting property of this class of models is that the energy band 
spectrum does not depend on a (see Fig. 1) whereas the eigenfunctions do. More quantitatively, 
the energy spectrum consists of two dispersive bands e±^k = a flat band at 

energy cq = A. For A > 0, the corresponding Bloch eigenfunctions (s = ±,0) read 
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Figure 1. (Color online). a-Ts model: (a) The dice or Ta lattice is constituted by two 
interpenetrating honeycomb lattices. Thick links, nearest-neighbor hoppings Cat from sites A 
to sites B {ca = Thin links, nearest-neighbor hoppings Sat from sites C to sites B 

{sa = By varying a, the model interpolates between honeycomb {a = 0) and dice 

(« = !). Potential term -|-A on sites A^C and —A on sites B. (b) Energy bands dispersion in 
k space, for A = 0: two dispersive bands and one flat band. This spectrum does not depend on 
a. 


where Ck = 


and Sk = 


l/fcl 


\/A2 + |/fe|2 + 

~ of each band is then readily obtained as 


with c| -|- s| = 1. The Berry connection 


Al = - 


1 — 

1 -|- 


Vfc0; 


^t = - 


1 + Cfc 




•^k — 


2(1 + Cfc) 




ki 


( 4 ) 


such that quite generally -|- A'^f^ + A~j^ = 0. From this Berry connection it is possible 
to obtain, for each band, the so-called Berry curvature x A^f^. This Berry curvature 

usually presents some strongly peaked structure near specihc k* points where interband coupling 
is strong, even if there is a gap separating the bands. Moreover, by considering a k-space 
closed orbit around such a k* point, it is possible to calculate a Berry phase that measures 

the winding of the phase 6k around such points k*. More precisely, considering closed orbits 
(e) that corresponds to a constant energy e = A^ -|- | /fc p, we deduce that for each band 
the Berry phase <1>|* = dk ■ Af, accumulated along such an orbit is given respectively by 

with = -4(1± A)$o^ where and dk-Vk6/2-K is 

the winding number of the orbit such that W<^ = 0 for a closed orbit around a regular fc*-point 
and = ±1 for an orbit encircling a Dirac point k* [7, 10]. 

From this perspective, the main interest of the a-Ts models described by Eq.(2) is that all 
these Berry quantities (connections, curvature and phase) are proportional to <1 >b = tt (see 
Eq.(4)). As a consequence, the interband effects that are encoded by these Berry quantities are 
reduced upon increasing the value of a from 0 and totally vanish for a = 1. Anticipating on 
what follows, we point out that our results for the orbital susceptibility of a-75 models provide 
clear evidence that all interband effects are not only encoded by the Berry quantities. 

























3. Orbital susceptibility formula 

For time reversal systems which we consider in the following, the orbital magnetization vanishes 
and the orbital susceptibility Xorh{P‘-,T) is the first measure of the sensitivity of the energy 
spectrum when it is placed in a uniform static perpendicular magnetic field B = Buz'. 
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where S is the area of the system, T is the temperature, /x is the chemical potential, ho = 
47r.lO“^S.I. and fl(r, |U, B) is the grand canonical potential of the non-interacting electrons gas 


Q{T,h,B) = -T 
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with p{u},B) the field dependent density of states (Dos) 

/9(w, 5) =--^mXr G(a; + i0+,5) (7) 

TT 

where G(a;, B) is the field dependent retarded Green’s function (we use /cb = 1 and h = 1). 

In order to obtain the susceptibility, it is thus necessary to calculate the density of states 
p{uj,B) or the Green’s function G{uj,B) = (cal — h{B))~^ associated to the field dependent 
Hamiltonian h{B). For the tight-binding model (1), h{B) is obtained by multiplying the real 
space hopping amplitudes by a gauge and position dependent Peierls phase: 


h{B) — '^^^ [Cat{6rg+5i,rAd-^rB+S2,rAd-^rB+S3,rA)^ b) {f" a\ /g\ 

~\~Scyt(^5rB—5\,rQ 5rB—S2,Tc — |^_B)(^d] “1“ h.C, 

where ‘Pry — e Ir ^ ' dZ with A{r) the gauge dependent vector potential associated to the 
uniform magnetic field. Starting from the Hamiltonian h{B) there are essentially two approaches 
to calculate p{uj,B). 

The first approach consists in computing the exact magnetic field dependent energy spectrum 
en{B) associated to h{B); this leads to the so-called Hofstadter butterfly spectrum. Using such 
an approach, we have recently studied the orbital susceptibility of the a-Ts model (2) for the case 
A = 0 [7]. In particular we have shown that the orbital susceptibility Xorh{H,T) strongly varies 
with the parameter a. More precisely, at the Dirac point, Xorhip- = 0,r) exhibits a continuous 
transition from a diamagnetic peak for a = 0 (honeycomb-graphene) to a paramagnetic peak for 
a = 1 (dice) (see Fig. 2). Away from the Dirac point, an opposite transition from paramagnetism 
to diamagnetism takes place such that a sum rule f d/xXorb(z^, T) = 0 is preserved for all a [7]. 
Despite these interesting results, such an approach suffers from being essentially numerical and 
thus it does not allow to fully understand how the dependence on the parameter a enters in the 
susceptibility. 

In order to better understand to role of the parameter a, we have developped a second 
approach which consists in calculating the susceptibility Xorb(zi) T) using our recently established 
gauge-invariant perturbative response formula [1]: 

Xorb = nF(u;) Tr{gh^-ghyy - gh^ygh^y - d{gh^gh^ghyghy - gh^ghygh^ghy)} , 

(9) 

where nF(w) = l/(e t -|- 1) is the Fermi function, g{id) = (cul — h)~^ is the retarded Green’s 
function associated to the zero-field Hamiltonian and = [x, Zi], = [x, [x, /i]] are the single 




Figure 2. (Color online). Orbital susceptibility obtained from the numerically 

computed Hofstadter spectrum [7]. x (iii units of the Landau band edge value \xl\ = 
(/ro/167r)(e^to^/^)) as a function of the chemical potential ^ (in units of t) in the whole band 
for various a as indicated and for a temperature T = 0.02t. At /i = 0 there is a transition from 
a diamagnetic peak for a = 0 (red: graphene) to a paramagnetic peak for a = 1 (blue: dice). 
Because of a sum rule, the orbital response at zero doping is systematically compensated by an 
opposite response at finite doping. 


and double commutators of the position operator with the zero-field Hamiltonian. As discussed 
in [1], the expression (9) is equivalent to a recent formula derived for graphene [11] but it differs 
from the well-known Fukuyama formula [6]. We stress however that when ^ 0, only formula 
(9) fully agrees with the susceptibility results obtained from direct numerical computation of the 
corresponding Hofstadter butterfly spectrum. Moreover, we remind that formula (9) is valid not 
only for Bloch electrons in infinite crystals but it also applies to disordered and finite systems 
as well. 

In the present paper, we restrict to multiband Bloch electrons in infinite crystals. In that 
situation, the trace operator is rewritten Tr(») = ti'(*) = S f ^|tr(») where the integration 

is performed over the first Brillouin zone (BZ) and tr(») is the partial trace operator on the 
band index. Accordingly, the susceptibility formula follows as 

f d'^k 

Xorh = — J J ^ nF{u})[Uk{u}) - 4I4(w)], (10) 


where 

Uk{oj) = ii{{gh-^ghyy - gh-ygh-y)k }, 

Vk{u:) = tv{{gh-gh-ghyghy - gh^ghygh^ghy)k} , ^ > 

with gk{^) = (wl — hk)~^ the Green’s function matrix associated to the zero-field Bloch 
Hamiltonian matrix and /i-^ = with (z, j) G (x, y). As detailed in the appendix, 

for the class of models described by Eq.(2), by defining the two components vectors (x = (xi, X 2 )) 
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Figure 3. (color online) Honeycomb lattice model that determines the form f^: t (full lines) 
and fs (dashed lines) respectively first and third nearest-neighbor hopping amplitudes. The two 
sites in the unit cell are called A and B and are shown as blue and red dots. In a uniaxially 
compressed honeycomb lattice, there are two different values for the nearest-neighbor hopping 
amplitudes: t for thin (non-vertical) lines and t' >t for thick (vertical) lines. 


it is possible to obtain the following compact expressions for Uk{uj) and Vkiuj): 


Uk{uj) = 2g+g-[{ul^uf 
Vk{uj) = {u%vl-ulvlf 


glgl[l+2>{\^)"]+Aglgl\fk\^ . 






(14) 

where ^±( 0 ;) = l/(a; — e±,k)- Interestingly, these last expressions show that Uk{oj) and Vk{oj) 
have poles only on the two dispersive bands as if the flat band did not play any role. In fact, 
the implicit effect of the flat band is essentially encoded in the a dependent term appearing in 
Vk{uj). However it is quite remarkable that Uk{uj) is totally independent of a. 


4. Orbital susceptibility of low-energy a-Ts models 

In this part of the paper we use equations (13,14) to calculate the susceptibility of a-Ts model for 
three different fk- In order to allow for a full analytical calculation, for each example we consider 
the low-energy effective model Hamiltonian that correctly describes the energy spectrum near 
the minimum of the corresponding \fk\- The three distinct forms of fk that we consider are 
defined on the honeycomb lattice model illustrated on Fig. 3, that corresponds to the a = 0 
limit in which the C sites are completely decoupled from A, B sites. The C sites are not shown 
on Fig. 3 for simplicity. 


4-1. a-Tz graphene model 

We first consider the usual isotropic model of graphene which corresponds to f = t and ts = 0 
in Fig. 3. In this situation, 

= t{e-^kya ^ cos{V3k^a/2)) (15) 

For A = 0, the corresponding a-Ta energy spectrum (see Fig. 1(b)) exhibits linear band touching 
at the two inequivalent corners ±K of the Brillouin zone. Introducing the valley index ^ = ±, 
the effective low-energy model is obtained by expanding fk near f^K+k — v{f,kx — iky) 

with the velocity v = ^ where a is the nearest-neighbor distance. Hereafter to simplify the 





notations of most equations we define the pseudo wavevector Kx,y = vkx,y For each valley, the 
low-energy a-Ts Hamiltonian reads: 
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Substituting these expressions in (10), summing over the two valleys and noting Kx 
K y = esin0, we obtain 
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Integrating over variables and following the procedure described in appendix B, we finally obtain 

1 — 
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which gives the correct behaviour for A = 0. For a = 0, Eq.( 21 ) with A = 0 recovers the 
famous diamagnetic peak originally found by McClure [ 12 ], whereas Eq.( 21 ) with A 7 ^ 0 
reproduces the in-gap diamagnetic plateau first obtained by Koshino and Ando [13]. Note 
however that these authors obtained their results by performing a low-field expansion of the 
grand potential calculated from direct summation over the effective Landau levels spectrum 

en{B) = ±eBY^[^ [12] or en{B) = A^ -|- e^|n| [13] where es = v\/2eB. For a / 0, 
Eq.(21) is also coherent with the results of such methods, albeit now with a-dependent effective 
Landau levels spectrum e„(i?) = ±£ 5 y^ln + 7 I or en{B) = A^ + e^ln + 7 I where n and 

7 = T+^ ['^]- picture, using the relation 7 = ^ — $B/27r where d>s is the a-dependent 

Berry phase of a Dirac cone [7, 10], the variation of the susceptibility with a is interpreted as 
a consequence of the variation of the corresponding Berry phase. As we already pointed out, 
for a = 1 the Berry phase vanishes and so do interband effects encoded by the Berry curvature. 
Nevertheless, for a = 1 there is still a paramagnetic susceptibility plateau inside the gap which 
is an indication that interband effects are still present. This result implies that some important 
interband effects are not encoded by the Berry curvature. More quantitatively we note that the 
a dependent prefactor of the susceptibility (21) changes sign at Oc = — 0.518. This signals 

a transition from a diamagnetic peak/p/ofeau for a < Oc to a paramagnetic peak/p/afeau for 
a > ccc- This value ac = 0.518 coincides with the numerical results obtained from computing 
the susceptibility with the Hofstadter spectrum [7]. As a final remark, we note that for the low- 
energy Hamiltonian (16) which is linear in h and thus separable (i.e. = 0), the Fukuyama 

formula [ 6 ] would have given the same results. 



4-2. a-Ts semi-Dirac model 

We now consider a model with t' = 2t and is = 0 (see Fig. 3). In that situation 


fk = cos{V3k^a/2)) 


( 22 ) 


For A = 0, the corresponding a-Tjj energy spectrum (see Fig. 4) exhibits a semi-Dirac band 
touching at one of the M points of the Brillouin zone [14]. The low-energy model obtained 

kP' 1 

by expanding near this point reads /ju-i-fe — ^ - ivky with the effective mass m* = ^ 

velocity v = ^ and it features a linear-quadratic spectrum Sk = =*= As in 

previous example, to simplify the notations of most equations we define the pseudo wavevectors 
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and Ky = vky. The low-energy a-Ts Hamiltonian is: 
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Quite interestingly, apart from some prefactor, the main noticable change between equations 
(20) and (27) is the effective density of states that appears in the integral: for the a-Ta semi- 
Dirac model it is oc whereas it is oc e in the a-Ts graphene model. From formulae (B.4,B.9) 
given in appendix B, we finally obtain: 

A = o,r = 0 

(28) 

Q+ |dl<A,T = 0. 

For 0 = 0 and A = 0, the first expression coincides with the result recently obtained in [1] 
for a two band model at the semi-Dirac point. For any a it predicts a diamagnetic peak that 
scales like 1 /bJfPnfJJxfT) in the gapless case A = 0. In the presence of a gap, the second 
expression shows that the susceptibility plateau in the gap scales as 1/VA. As in the previous 
example, Eq.(28) predicts a diamagnetic to paramagnetic transition for a > cxc at a critical value 
ac = \/2 — 1 — 0.414 which is smaller than in then previous example. We have verified that this 
value ac = 0.414 agrees with the numerical results obtained from computing the susceptibility 
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Figure 4. (Color online). a-Ts semi-Dirac model energy bands dispersion in k space, for A = 0. 
(a) Full spectrum with semi-Dirac band touching of the two dispersive bands, (b) low-energy 
spectrum. 



with the Hofstadter spectrum. We stress that for the a-Ts semi-Dirac model, as yet, there is no 
exact analytical expression for the Landau levels e„(i?) of the low-energy model. The Landau 
levels are related to the eigenvalues of a modified quartic oscillator [14]. However by computing 
the Hofstadter spectrum of this a-Tjj semi-Dirac model for small magnetic field and various a, 
we have observed that the effective Landau levels en{B) are well described by the approximate 
form €n oc [(n -|- 1/2)H]^/^, first obtained in [14]. Interestingly, for n > 3 the numerical Landau 
levels en{B) appear to be almost independent of a. This observation implies that the dia- to 
paramagnetic transition is essentially driven by the variation of the n = 0,1, 2 Landau levels. We 
also stress that for this semi-Dirac model, there is no Berry phase around the semi-Dirac point. 
However the Berry curvature is non-zero (for a 7^ 1) and it exhibits a double-peak structure 
near the semi-Dirac point. There is still a lack of a clear physical picture of the origin of the 
dia to paramagnetic transition for this a-Ts semi-Dirac model. As a final remark, we note that 
since the low-energy Hamiltonian (23) is still separable (i.e. = 0), the Fukuyama formula [6] 

would have given the same result. 


4-3. a-Tz pseudo-bilayer model 

The last model we consider has parameters t' = t and = 1/2 in Fig. 5 [15]. In that situation, 
/fc = cos(\/3fc^a/2)) + ^ + 26"*^^“ cos{V3k^a)) (29) 


For A = 0, the corresponding a-Ts energy spectrum (see Fig. 5) exhibits quadratic band touching 
at zizK similar to a bilayer graphene. The low-energy model obtained by expanding near 
reads f^K+k — ~ with the effective mass m* = ^ and the energy spectrum 


£k = 2 ^- in previous sections, to simplify the notations of most equations, we define the 
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The low-energy a-Ts Hamiltonian becomes: 
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Figure 5. (Color online). a-Ts pseudo-bilayer model energy bands dispersion in k space, for 
A = 0. (a) Full energy spectrum with quadratic band touching, (b) low-energy spectrum. 


This effective Hamiltonian is now quadratic in Kx^y and it is not separable (i.e. 7 ^ 0). As a 

consequence, we do not expect the Fukuyama formula to give the correct result in that situation. 
Following similar steps as in previous examples we obtain 


f={Kl-Kl,-2^KxKy), = 2{Kx,-CKy), P = 2{-Ky,-^Kx), 

p - = 2(1,0), fyy = 2(-i, 0), py = 2(0, p ), 


such that by noting Kx = ^/scos6, Ky = ^/ssm0 we deduce 
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Substituting these expressions in (10), summing over the two valley, integrating over 9 we find 
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From Eqs.(B.4,B.9) given in appendix B, we finally obtain: 
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where it is necessary to introduce a cutoff scale Cc [16, 17]. As in previous examples, 
there is a diamagnetic to paramagnetic transition of the logarithmic term at a critical value 
Oc = \/4 — \/3 = 0.263, whereas there is a paramagnetic to diamagnetic transition of the second 
term at ac = \/l2 — \/TT = 0.147. For a = 0, these expressions coincide with previous results 
obtained in [16, 17]. In [17], the susceptibility was derived from using the analytical form 
of the effective Landau levels = eB\/n{n — 1) (e^ = ^) that are associated to the low- 
energy model. A similar derivation for any a seems difficult to achieve even if the Landau level 























sequence is known. In that perspective, we note that for the gapless a-Ts pseudo-bilayer model 
it is possible to derive a topological Berry phase that varies with a; this already gives some 
insight on the modified semiclassical Landau levels spectrum [10]. However as noted in [10], 
for the bilayer system the semiclassical Landau levels spectrum does not fully agree with the 
exact quantum Landau levels spectrum and therefore we do not expect the semiclassical Landau 
levels spectrum to be sufficient to describe correctly the susceptibility -especially the singular 
behaviour of the susceptibility near half filling /x = 0. 

5. Summary 

In this work we have studied the importance of interband effects on the orbital susceptibility of 
three bands a-Ts tight-binding models that were initiated in [7]. The particularity of the a-Ta 
tight-binding models is that the coupling between the three energy bands (which is encoded 
in the wavefunctions properties) can be tuned (by a parameter a) without any modification of 
the energy spectrum. To highlight the role of these interbands effects, we have examined the 
orbital susceptibility of these models. Using the gauge-invariant perturbative formalism that we 
have recently developped [1], we obtained the generic formula of the orbital susceptibility of the 
a-Ta tight-binding models. More quantitatively we have calculated the orbital susceptibility of 
three distinct a-Ta tight-binding models: a-Ta graphene model, the a-Ta semi-Dirac model and 
the a-Ta pseudo-bilayer model. To obtain an analytical form of the susceptibility, we have only 
considered a low-energy Hamiltonian of these models; it correctly describes the energy spectrum 
near the band touching where it is expected that interband effects are the strongest. The main 
result of our study is that for each of these models, by varying the parameter a and thus the 
wavefunctions interband coupling, it is possible to drive a transition from a diamagnetic to a 
paramagnetic behavior, both in absence or presence of a gap separating the dispersive bands. 
In particular we emphasize that the in-gap susceptibility does not need to be diamagnetic. 
Moreover the existence of a finite in-gap (paramagnetic) susceptibility at a = 1 (for each model) 
provides hints that some important interband effects are not encoded in the Berry curvature 
which vanishes for a = 1. 


Appendix A. Determination of Uk{oj) and Vk{oj) 

The aim of this section is to give the key step to find the expressions Uk{u}) and of Eq.(14). 

We first slightly rewrite their definitions Eq.(ll): 


Uk{uj) = iT{{gh--ghyy - gh-ygh-y)k }, 

Vkiu;)=tT{{gh-ghy[ghy,gh^])k}, 


where gk{<-^) = (^^I — hk)~^ it the Green’s function matrix associated to the zero-field Bloch 
Hamiltonian matrix Eq.(2), h], = and with (i,j) = {x,y). The calculation of 

Uk: hfc is made simple by remarking that the matrices h, g, can each be written in terms 

of the following three matrices only: 




T 

-H 

o 

0 \ 


( 1 

0 

0 \ 

Co-e*®'® 0 


II 

0 

-1 

0 

\ 0 Sae^^'^ 

0 J 


V 0 

0 



(A.2) 


Note that the two matrices S± depend on both parameter a and wavevector k; moreover we 
stress that S- is an antihermitian matrix whereas are hermitian. With these definitions, it 
is immediate to verify the identities 


h = ASo + \fk\S+, 
h[= uiS+ + iv{S-, 
= Uk^+ -h wiS-, 


(A.3) 



with 


u 

u 


j _ fj-fk 
\fh\J 

k - l/fel ’ 



k 


ifkXpk) 

l/fel 

ifkXfl^) 
■ l/fel 


where we have defined the two components vectors {x = (xi,X 2 )) 


fk = {^efk,^rnfk), 

B _ f 9SRe/fc d^mfk n 
fc i 9fci ’ 0fci f ’ 
A’ _ /- Q^SRe/fc d^'imfk \ 
V dkidkj ’ dkidkj 


(A.4) 


(A.5) 


To obtain a simple form for g we remark that the matrices appear to have very peculiar 
properties: = ±5^, Sq = I, [5±,5o]+ = 0. From these properties we deduce that 

g = (cjI — hk)~^ can be written as g = ail + a 2 So + 03 ^+ + 04 ^^ where 01 , 2 , 3,4 are determined 
from using the identity (uj — h)g = 1. We then obtain: 

a = 5+5-M - ^Sq - |/fc|S'+) - g+g-go\fk\^{l - 5+), (A.6) 


where g±k{oj) = ^ and go{oj) = . Note that g{oj) has three poles corresponding to 

the three bands. Substituting identities Eqs.(A.3,A.6) in Eq.(A.l), we obtain the expressions 
Eq.(14): 


C/fc(a;) = 2g+g_[{ul^uf 
Vkiuj) = {ulvl-ulvlf 


- ~ + ^a+a-{uk^u^k 

5i<7^[l + 3(^)']+44ff^l/fcP ■ 






(A.7) 


Appendix B. Pole decomposition and integration over variables e and uj 

The generic form of the susceptibility of the three different a-Tjj models is: 


Xorb 


oc^fduj npiuj) de u(e) 
oc^fduj nF(uj) /q“ de u(£) 


Ae^ I , C 

(a;2_(A2+£2))3 -b (^2_(A2+£2))2 (^2_(A2+£2)) 

[Ae'^{g+g-f + B£‘^{g+g_f + C{g+g-)] 


(B.l) 


where gs{oj) = l/(w — e^) and with e* = + e^ {s = ±). The effective density of states 

12 (e) and the parameters A, B, C of the three models (a-graphene, a-semi-Dirac, a-bilayer) are 
summarized in table Bl. 



a-graphene 

a-semi-Dirac 

a-bilayer 

12(e) 

11 £ 

l/^/^ 

1 

A 

4 

4 

16 

B 

i+Kiai) 

l + 3 (^) 

5 + 12 ( 1 +^) 

C 

0 

0 

1 


Table Bl. Effective density of states iu{£) and parameters A,B,C of the three considered 
models: a-graphene, a-semi-Dirac and a-bilayer. 



We now describe the key steps to perform the explicit integration over variables e and oj. 
The first step consists to separate the poles at Eg'- 

Xorb oc ^ J du npiuj) de u(e) + (f fj “ + (^ “ f If + 

(B.2)' 

The next step consists to use the identity ^ f dcj nF{uj)g^{uj) = — where 

is the derivative; we can the rewrite: 


Xorb OC — X]s=± fo + (f fj “ 3^||)np(es) + — f 

(B.3) 

From this point, we separate the study of case A = 0 from case A 7 ^ 0. 

Gapless models, A = 0; 

For A = 0, since = se we first rewrite: 

Xorb oc - /o“ de iy(£) ^ein'pie) - np(-e)) + (f - 34)(f^F(e) + np(-e)) + (§ - f + 3 ^)!AM:^At£) 

(B.4) 

We now note that for the three considered cases, i^(e) ~ e^ with — 1 < p < 0; this property 
permits an integration by part of the term proportionnal to Up and np: 


A / 


/q“ deeP+i(rap(e) - np(-e)) = 
deeP"i(nF(e) - nF(-e)) = 


/ / 00 . . 

£P''~^(np(£) + np(-e)) - (p + 1) deeP(np(e) + np(-e)). 


£P 


,oJ0 


p (nF(e)-nF(-e))J^ -/“ de^(np(e) + np(-e)). 


(B.5) 

where the last line requires p < 0. Using these identities for the three considered cases 
(p = —1, —1/2, 0) and using the parameters of table B1 we deduce: 


Xorb oc - 5]np(0), 

,, 3r/ 1 -| poo ^^np(s)+n‘ 

Xorb OC 2 JJ 0 vT 


Xorb OC 3 


1(1 


_ 3j joo nF(£)-nF(-£) ff 1-a^ ^ _ 11 


l+a^J 


+ 


l\l+o2 J I 2 J 


a—graphene 
a—semi — Dirac 
a—bilayer 


(B.6) 


Gapped models, A / 0; 


For the gapped case, we can also perform some integration by part of terms proportionnal to 


drEl 


rip and np. More precisely using the identities and ^ ^ we first obtain 


rd£^n; =-/“de^(p + 3-2|J; 


rOO ^ ^P+2 

fo deV 


F 

np 


n 




pp+i 

— —np 


TOO 1 pP+4 

Jo 


gP+3 


-lo. 


-np 


+ r de 
'^+ird£ 


{p+l)eP 


np + 


fP+2 


-nr 


(p+3)£ 


P+2 


2 L 

pP+3 

_ Tl 7-1 

Jo 

OO , „ 

1 P+3 

^ L 

pP+i 

_ ^ _ Tl 7-1 


0 



OO 


-np + 


+ Uo^d£ 


_p+4 / 

{p+3)(p+l)sP 


I {p+3)ePe^ 

np + -n 


rP+4 


F T 

(B.7) 


-n. 



from which we deduce: 


Xorb OC X]s=± 


' A 

£P+3 

16 

-^np 


J 0 


+ [f “ 


ep+1 


-np 


Jo 


+(f - (P + l)[f - t|(p + 3)])/o“ de fjnF(es) 


(B.8) 


For the three considered cases {p = —1,—1/2,0) with the parameters of table B1 we finally 
obtain; 


Xorb OC 
Xorb OC 
Xorb OC 


3r/l-a2\2 

4 l+a '-^ J 


l-inp{A)-np{-A) 
3 J A 


3r/_ 1] foo i^ np(\/A^+£^)-np(-\/A^+e^) 
8\-\T+^ J 2\Ja v^VA2+£2 

/ 1-q 2'^2 j^np(VA2+e)-nF(-VA2+£) 

J 4JJa V^A^ 




cr—graphene, 
a—semi — Dirac 
a—bilayer 

(B.9) 
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